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It is shown that Majorana fermions trapped in three vortices in a p-wave superfluid form a qubit 
in a topological quantum computing (TQC). Several similar ideas have already been proposed: 
Ivanov [Phys. Rev. Lett. 86, 268 (2001)] and Zhang et al. [Phys. Rev. Lett. 99, 220502 (2007)] 
have proposed schemes in which a qubit is implemented with two and four Majorana fermions, 
respectively, where a qubit operation is performed by exchanging the positions of Majorana fermions. 
The set of gates thus obtained is a discrete subset of the relevant unitary group. We propose, in this 
paper, a new scheme, where three Majorana fermions form a qubit. We show that continuous 1- 
qubit gate operations are possible by exchanging the positions of Majorana fermions complemented 
with dynamical phase change. 2-qubit gates are realized through the use of the coupling between 
Majorana fermions of different qubits. 
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-f— > 



> 

\6 
o 
o> 
o 



X 



I. INTRODUCTION 

Ivanov first pointed out that a pair of Majorana 
fermions can be used to implement a qubit and proposed 
gate operations on it [l[ . He has also demonstrated that 
a braiding of Majorana fermions leads to entanglement 
of two qubits. Later, Zhang et al. proposed to use four 
Majorana fermions to implement a qubit Q. They fur- 
ther proposed to use a flying qubit to entangle two qubits 
thus implemented. It should be noted, however, that a 
braiding is a discrete operation and it is impossible to 
implement an arbitrary one-qubit gate with a braiding. 
Moreover, it should be also pointed out that entangling 
operation using a flying qubit does not work in practice, 
since the Majorana fermion does not couple with density 
fluctuation as shown in Q. It is the purpose of this pa- 
per to show that continuous gate operations are possible 
if a qubit is implemented with three Majorana fermions. 
We use two Majorana fermions, similarly to Ivanov's pro- 
posal, to implement a qubit and an additional Majorana 
fermion for continuous control of the qubit state. Sim- 
ilarly continuous 2-qubit gates can be implemented by 
making use of the coupling between Majorana fermions 
which belong to different qubits. 

Let us consider a p-wave superfluid with the order pa- 
rameter p x + ipy A vortex in the superfluid supports a 
bound state in the quasiparticle spectrum, whose bound 
state energy is exactly at the center of the band gap. 
The bound state is invariant under charge conjugation 
and called the Majorana mode, which will be called the 
Majorana fermion hereafter [4]. It has been shown by 
Mizushima, Ichioka and Machida that this zero-energy 
state is energetically well separated from the other bound 



states (Caroli-de Gennes-Matericon states) in the strong 
coupling limit, in which the energy gap A is on the same 
order as the Fermi energy Ep |5|. Topological quantum 
computing employs Majorana fermions in such strongly 
coupled systems p. 

Let us consider a two-Majorana fermion system, first. 
The Hamiltonian of this system is given by 



H = i J127172, 



(1) 



where J12 is the coupling constant between two Majo- 
rana fermions and ji stands for the Majorana operator 
associated with the ith vortex. They satisfy the anticom- 
mutation relation 

{ 7i , 7j } = 2S ij . (2) 
We now introduce another set of operators a and 

a = -(71 + 172), a f = -(71 - 172), (3) 

which satisfy the fermion anticommutation relation 

{a, a} = {a\cJ} = 0, {a, a f } = 1. (4) 

The Hamiltonian is then rewritten, in terms of the new 
operatos, as 



H = w (2a t Q ! - 1) 



(5) 
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It is shown that the Bogoliubov wave functions u(r) 
and v(r) satisfy the relation u(r) = v*(r) for a zero- 
energy mode and hence the Majorana operator is ex- 
pressed as 7^ = Q + cj, where q = J d 2 rui(r)*ip(r). Here 
ip(r) is the field operator of the particles in p-wave su- 
perfluid state and Ui(r) is the Bogoliubov wave function 
of the zero-energy state trapped in the ith vortex. Let 
|0)i, defined by Cj|0)j = 0, denote the state in which the 
ith vortex has no zero-energy particle, while c]|0)j = |1); 
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denote the state with a zero-energy particle at the ith 
vortex. 

The ground state energy of this Hamiltonian ([5]) is — u), 
which has two eigenvectors 



a|0)i|0) 2 , aa t |0)i|0) 2 , 



(6) 



where the first eigenvector has odd fermion number while 
the second one has even fermion number. The excited 
state energy is +oj, which is doubly degenerate with the 
energy eigenstates 



a T |0)i|0) 2 , ata|0)i|0) 



(7) 



where the first eigenvector again has odd fermion number 
while the second one has even fermion number. 

We note that application of on the ground states 
changes the parity of the fermion number as 



a|0)i|0) 2 -> a ta|0)i|0) 2 
aat|0) x |0)2 -> a t |0) 1 |0) 2 . 



(8) 



II. THREE-MAJORANA FERMION MODEL 

Suppose there are three vortices, each of which sup- 
ports a Majorana fermion in the limit of infinite sepa- 
rations among the vortices. The Hamiltonian describing 
the coupled Majorana fermions is given by 



H = lJl2lll2 + ^237273 + ^317371; 



(9) 



where J„ £ R is the coupling strength between the ith 
and the jth Majorana fermions. It turns out to be con- 
venient to parametrize three coupling constants by the 
polar angles 6 and as 

(■^23) J31, J12) — J {sin 8 cos (j), sin 9 sin (j), cos 9). (10) 

The Hamiltonian is diagonalized by introducing the 
creation and the destruction operators 



vt — I 



[(cos 9 cos cf) + i sin 0)71 

+ (cos 9 sin <fi — i cos 4>)l2 

ai=\ [(cos 9 cos — i sin 0)71 

+ (cos 9 sin 4> + i cos (j))^ 



sin 6*73] 
sin 073] 



(11) 



as 



H = u (2. 



1 



(12) 



where oj = J. It should be noted that there exists a 
Majorana fermion operator 



ft = sin 9 cos <pji + sin 9 sin <f>j2 + cos $73 , 



(13) 



which is orthogonal to a and cr . It is easy to verify these 
fermionic operators satisfy the anticommutation relations 



It follows from the above anticommutation relations that 
ft commutes with H and, hence, ft represents the zero- 
energy Majorana fermion. Mizushima and Machida 
analyzed the lowest energy eigenvalues by solving the 
Bogoliubov-de Gennes equation numerically and ob- 
tained the same results @- 

The operators a, or and ft take the simpler forms 



-i<j> 



i<l> 



a = ^—-(71 +i7 2 ),a t = ^-(71 - £72), £ = 73 (15) 



in the limit Ji 2 3> J 2 3, J31, which corresponds to the case 
in which vortex 3 is isolated from vortices 1 and 2. We 
also have 



tan m = 



J31 

</ 2 3 



(16) 



Now let us analyze the energy eigenstates of the Hamil- 
tonian © in the above limit. The ground state with the 
energy — u> is four-fold degenerate. Ground states with 
odd number of Majorana fermions are two-fold degener- 
ate, 



a|0)i|0) 2 |0) 3 = 



-i(f> 



-i<f> 



-(7i+£72)|0)i|0) 2 |0) 3 
(|1) 1 |0) 2 |0) 3 + £|0) 1 |1) 2 |0) 3 ) 



(17) 



aa t /3|0) 1 |0) 2 |0) 3 



= ~ *7i7273)|0)i|0) 2 |0) 3 

= § (|0) 1 |0) 2 |1) 3 - i!l)ill) 2 |l) 3 ) • 



(18) 



Similarly, ground states with even number of Majorana 
fermions are two-fold degenerate with the eigenstates 



aa t |0) 1 |0) 2 |0) 3 = 



-(l-£7i72)|0)i|0) 2 |0) 3 
^(|0)i|0) 2 |0) 3 -£|l) 1 |l) 2 |0) 3 ) 



(19) 



-i<p 



*/3|0)i|0) 2 |0) 3 = — (7i+*72)73|0)i|0) 2 |0) 3 



-i<f> 



(|l)i|0) 2 |l) 3 + i|0)i|l) 2 |l) 3 ). 



(20) 



{a,at} = l, {a,ft} = {a\ft} = 0. 



(14) 



The excited state with the energy w is also four-fold 
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degenerate; states with odd fermion number are 



III. ONE-QUBIT GATES 



a t |0)i|0) 2 |0) 3 = ^(71 -*Ta)|0)i|0) 3 |0) 3 



i<l> 



Let us first consider the odd fermion number sector 
with the initial state 



(|l)i|0) a |0) 3 -t|0)i|l) a |0) 8 ) 



(21) 



a|0) 1 |0) 2 |0) 3 



-i(f> 



1 



-i(f) 



a T a/3|0) 1 |0) 2 |0) 3 = -(73 + «7i727 3 )|0)i|0) 2 |0) 3 



-( 7l +z 72 )|0)i|0) 2 |0) 3 
-(|l)i|0) 2 |0) 3 + z|0) 1 |f) 2 |0) 3 ). 



(|0)i|0) 2 |l) 8 + i|l)i|l) 2 |l) 3 ) 



(22) 



while those with even fermion numbers are 
a t a|0) 1 |0) 2 |0) 3 = |(1 + «TiT2)|0)i|0) 2 |0) 3 

= 5(|o)i|o) a |o) 3 + i|i)i|i) a |o) 3 ) 



(23) 



t /3|0) 1 |0) 2 |0) 3 = — (7i-n2)73|Q)i|0) 2 |0) 3 



■(|l)i|0) 2 |l) 3 -i|0)i|l) 2 |l) 3 ). 



We assume the vortices at 1 and 2 are also remotely sepa- 
rated initially so that all the coupling strengths are small. 
We still impose the condition J12 3> J23, J31 even in this 
case. Then the dynamical phase changes for the ground 
states and the excited states are almost identical since 
lj = J is negligibly small. Now we outline how to im- 
plement a unitary gate with continuous parameters in 
several steps as shown in Fig. 1. 

STEP f Suppose vortices at positions 3 and 1 are ex- 
changed in the counterclockwise sense, as shown 
in Fig. 1 (a), so that the Majorana operators are 
transformed as 73 — > 71 and 71 — > — 73. Under this 
transformation, the operator a transforms as 



(24) 

Transitions among the ground states and the excited 
states could be performed by Rabi oscillation through 
modulation in J 23 or J 31 . Suppose the interactions 

iSJ 2S j2jsCOs2ujt = -5J 23 (o!+e-^ - ae^)/3cos2w£ 

i5J 3 i7 3 7i cos2wi = — iSJ3i(a^e~ l< ^ + ae l ^)/3 cos 2oot 

(25) 

are introduced in the Hamiltonian. Then the follow- 
ing Rabi oscillations take place between the four sets of 
states; 

ground state excited state 

a|0)i|0) a |0) 3 o ata£|0)i|0) 2 |0) 3 
aat/3|0)i|0) 2 |0) 3 o at|0)i|0) 2 |0) 3 (26 ) 
«« t |0)i|0) 2 |0) 3 o atj3|0) 1 |0) 2 |0) a 
^|0)i|0) 2 |0) 3 o at a | ) 1 |0) 2 |0) 3 . 

Note that the Rabi oscillations preserve the parity of the 
fermion number. It is possible to implement a continuous 
series of quantum gate operations by making use of the 
above Rabi oscillations. However, this may cause qubit 
operation error since the system is under external field, 
which possibly contains noise. It is certainly desirable to 
perform qubit operations without errors by exchanging 
the vortex positions as was proposed by Ivanov |l| and 
Zhang et al. 0. 

Now we turn to our main result, in which continuous 
qubit operations are implemented by introducing dynam- 
ical phases in TQC. 



-i(f> 



-(71 + ij 2 ) 



-(-73 + «72) 



-i<j> 



-(aa f /3 + a ] af3) 



-i<j> 



-(ae 



i<j> 



) 



Transformations of the operators a\aa*(3 and 
a) uf3 under this exchange are also obtained and 
summarized as 



(27) 
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where 



m 3 i 



/ 1 

e i<t> 

e i<t> 



\ 



a 2i<j> 



,i<t> 



-i(f> 



ict> 



-2i<j> 
-i(j> 



(28) 
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STEP 2 Vortices at I and 2 are put close to each other, as 
shown in Fig. 1 (b), so that Ji 2 is appreciably large. 
Now both the ground state and the excited states 
acquire nontrivial phases. The transformation ma- 
trix is 



/ e- ir > 
e ir > 


V 



\ 





e-"> 




(29) 



/ 
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(a) 



• 1 

(b) 




(c) 



FIG. 1: Implementation of a one-qubit gate. Numbers 1, 2 and 3 show the positions of vortices, (a) Vortices at positions 1 
and 3 are exchanged (STEP 1). (b) Vortices at 1 and 2 are put close to each other so that they acquire the dynamical phase 
(STEP 2). (c) Vortices at 1 and 3 are exchanged again so that the vortices take their initial configuration (STEP 3). 



STEP 3 Subsequently, vortices at 3 and 1 are exchanged 
in clockwise sense as shown in Fig. 1 (c), which 
introduces rn^. 

The above three steps result in a transformation matrix 



m 31 m z m 31 
( 



COST] 


— ie ^ sin 77 








ie %< ^ sin 77 


cos 77 














cos 77 


ie~ lc ^ sin 77 








ie 1 ^ sin 77 


cos 77 



\ 

/ 

(30) 



This result shows that the qubit basis vectors |0) = 
a|0)i|0) 2 |0) 3 and |1) = a t a/3|0)i|0) 2 |0) 3 are continuously 
transformed. This statement remains true if another set 
of the qubit basis vectors, |0) = aa 1 '/3|0)i|0)2|0)3 and 
|1> = at|0)i|0) 2 |0) 3 , are chosen. 
It is instructive to implement the Hadamard gate 



1 



1 1 



V2 V 1 - 1 



with our scheme. We use |0) = o:|0) 1 1 0)2 10)3 and |1) = 
q;^q;/3|0)i |0)2|0)3 as the qubit basis. Then the upper-left 
block of the matrix (|3T))) has relevance. Let us write 



M(r,,ct>) 



COST] 

-ie 1 '*' sin 77 



-ie ^ sin 77 
cos 77 



(31) 



Then we easily verify the product M(f , -f )M(§, 0) im- 
plements the Hadamard gate up to an overall phase. 

Qubit operations are also possible by exchanging vor- 
tices at 2 and 3, instead of vortices at 1 and 2. It is also 
easy to verify that a similar qubit construction and qubit 
operations are possible if the qubit basis states are made 
of even fermion number states. The sequence of opera- 
tions given in Fig. 1, in this case, results in the matrix 
(f3"0")) . although m 3 i takes a different form from the odd 
fermion case (l28l) . 

It has been shown so far that a continuous family of 
1-qubit operations can be implemented by adding a third 
Majorana fermion to a pair of Majorana fermions. 



IV. TWO-QUBIT GATES 

Finally, we show that our qubits satisfy the universal- 
ity criterion by demonstrating that two-qubit gates can 
be implemented within the current proposal. We first 
note that the third Majorana fermion is required only to 
implement single-qubit gates and plays no role if it is far 
remote from the first and the second Majorana fermions. 
Let us first consider the braiding proposed in 1]. Let 71 
and 72 ('y'l and 7 2 ) be the Majorana fermion operators 
associated with qubit 1 (2), where an index associated 
with the second qubit is denoted with a prime. Let the 
initial state of qubits 1 and 2 be aa?\0), where 



a = -(7! +772): 



and we write |0) 1 1 0)2 1 0) i' , 1 0)2' as |0) to simplify the nota- 
tion. Ivanov attempted to create an entangled state 
^7|(|0)|0) + by braiding of Majorana fermions. Let 
us exchange Majorana fermions 1 and 1' in the counter- 
clockwise sense. The state then transforms as 

aa'\0) ->• l;{aa' + a^a'^ - aa'a'a 4 + aJaa'^a')\0), 

which is certainly an entangled state. However, this state 
is different from the state 

i(|0>|0) + |1>|1» = l(oa' + aWV)|0), (32) 

for example, to be implemented 

Now we would like to propose an alternative opera- 
tion to implement the state (|3"2"|) . We first let Majorana 
fermion 71 of qubit 1 and Majorana fermion "f[ of qubit 
2 come closer so that they interact with each other. The 
relevant interaction Hamiltonian is 

iJivJiJv — iJ\v (a + or) (a' + a''). (33) 

The interaction strengths are arranged to satisfy 

|Ji 2 |,|Ji'2'|»|Jii'| (34) 
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FIG. 2: Two-qubit system. Physical quantities associated 
with the second qubit are denoted with a prime. The coupling 
strength between Majorana fermions 1 and 1' is denoted as 
Jxi', for example. 

and 

|Jii'|»|Jia-^i'2'|. (35) 

It follows from the condition that the state aa'|0) 
has no time evolution since J\y is negligible compared 
to J\2 + Ji'2' • In contrast, there is an oscillation between 
two states aa'^a'\0) and a^aa'\0) since it follows from 
the condition (|35p that | J\2~ Jl'2' I is negligible compared 
to Jul Now we are ready to outline how to generate a 
state like (l33l) . 

STEP 1 We first prepare the state aa'\0). 

STEP 2 Apply M(tt/2, tt/2) of Eq. ([31]) on the second qubit 
to generate a state (W'a'lO). 

STEP 3 Introduce J\y coupling to transform the state into 

-= (aa n a' +aW) 10). 
y/2 

STEP 4 Apply M(tt/2,tt/2) again on the second qubit to 
obtain the entangled state 

—/={— olol + aftaafta')\Q) (36) 



as promised. 

We have dropped the operators and /3' which appear 
in the intermediate state. 

There is practically no change in the state (|36|) due to 
the condition once this state is created. Qubits 1 
and 2 may be widely separated for further stabilization. 



V. CONCLUSION 

In conclusion, we have proposed new qubit construc- 
tion in topological quantum computing, in which Majo- 
rana fermions trapped in a two-dimensional p-wave su- 
pcrfluid are employed. A single qubit is constructed out 
of three Majorana fermions. An arbitrary one-qubit gate 
can be implemented by a combination of the braiding of 
the vortices (and hence the Majorana fermions) and the 
dynamical phase change. Entangling operation required 
for two-qubit gate implementation is shown be realizable 
in a similar manner. 

Introducing a dynamical phase in TQC might seem to 
be a flaw in an otherwise perfect quantum computation 
scheme. It should be noted, however, that a brading in 
mathematics, which requires exact exchange of positions 
of Majorana fermions, is never possible to realize physi- 
cally. Exchange of positions in reality always involve an 
imperfection. 
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